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Measurement of photon-number statistics of fields com- 
posed of photon pairs generated in spontaneous paramet- 
ric downcoversion pumped by strong ultrashort pulses is de- 
scribed. Final detection quantum efficiencies, noises as well 
as possible loss of one or both photons from a pair are taken 
into account. Measured data provided by an intensified single- 
photon CCD camera are analyzed along the developed model. 
The joint signal-idler photon-number distribution is obtained 
using the expectation maximization algorithm. Covariance of 
the signal and idler photon-numbers equals 80 %. Statistics 
of the generated photon pairs are identified to be Poissonian 
in our case. Distribution of the integrated intensities of the 
signal and idler fields shows strong correlations between the 
fields. Negative values of this distribution occurring in some 
regions clearly demonstrate a nonclassical character of the 
light composed of photon pairs. 

I. INTRODUCTION 

Light generated in the process of spontaneous paramet- 
ric downconversion is emitted in photon pairs [1]. Pho- 
tons comprising one photon pair are strongly correlated 
— they are entangled. Entanglement of photons in a pair 
has been used in many experiments that have provided 
a deep insight into the laws of quantum mechanics [2,3]. 
Among others, the measured violation of Bell inequali- 
ties ruled out neoclassical local hidden- variables theories. 
Photon pairs have also found ways to practical appli- 
cations, e.g., in quantum cryptography, measurement of 
ultrashort time delays, or absolute measurements of de- 
tection quantum efficiencies. These experiments utilized 
photon fields that contain only one photon pair in a mea- 
sured time window with high probability. 

There have been experiments (teleportation, measure- 
ment of GHZ correlations, etc.) measuring triple and 
quadruple coincidence counts caused by fields contain- 
ing two photon pairs in the time window given by an 
ultrashort pump pulse during the last couple of years. 
However, states used for such experiments contain a very 
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low fraction of two photon-pair states in comparison with 
one photon-pair state and the vacuum state. The rea- 
son is to eliminate the influence of three and more-than- 
three photon-pair states in considered experimental se- 
tups. Measurements done in such setups are conditional 
and they require long data-acquisition times. 

The use of more powerful pulses pump lasers as well as 
development of materials with higher values of x' 2 ' sus- 
ceptibilities open the way to generate fields containing 
more than one photon pair from one pump pulse with 
nonnegligible probabilities. For such fields, the influence 
of states containing two photon-pairs has to be judged 
even in experiments utilizing primarily just one photon 
pair. For example, the effect of two-photon-pair states 
in measurement of time-bin entanglement has been an- 
alyzed in [4]. Contribution from more-than-one photon- 
pair states becomes essential for higher mean values of 
photon pairs and such fields then have to be character- 
ized in general by their photon-pair statistics. 

Determination of photon-pair statistics is necessary 
also for characterization of weak continuous downcon- 
verted fields if they are detected in long time windows 
[5]. In this case photon-pair statistics were determined 
to be Poissonian if dead-time effects in their detection 
were eliminated [5] . 

We describe the behaviour of photon pairs in a general 
photon-number resolving system. Assuming detection of 
generated photon pairs by an intensified CCD camera, 
we can obtain a joint photon-number distribution of the 
signal and idler fields at the output plane of a crystal 
from measured data using a developed model of detection 
and a suitable reconstruction algorithm. The generated 
field can then be characterized in general by the joint 
distribution of the signal and idler integrated intensities. 

The paper is organized as follows. Sec. 2 contains 
a general model describing a photon-number-resolving 
detection device. Sec. 3 is devoted to reconstruc- 
tion of the joint signal-idler photon-number distribution. 
Joint signal-idler integrated-intensity distribution is de- 
termined in Sec. 4. Sec. 5 provides conclusions. Statis- 
tical properties of fields generated in spontaneous para- 
metric downconversion are derived in Appendix A from 
first principles. 
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II. MODEL DESCRIBING THE MEASUREMENT 
OF MULTIPHOTON-COINCIDENCE COUNTS 

Measurement of the joint signal-idler photon-number 
distribution can be in general described using the scheme 
shown in Fig. 1 [6-8]. 



pump field 




1 x N s multiport 



FIG. 1. Scheme of the considered model. Photon pairs 
are generated in nonlinear crystal NLC. Virtual beamsplitters 
BSs and BSj describe possible losses of one or both photons 
from a pair before they are detected. Signal (idler) photons 
propagate through a 1 x Ns (1 x iVj) multiport and are de- 
tected in single-photon detectors Ds 1; Ds 2 , . . . , £>s Ns (D/ 1; 
D/ 2 , . . . , T)i N ). Signals from the detectors are registered in 
a coincidence-detection device. 

Photon pairs occurring in the output plane of nonlinear 
crystal NLC propagate towards photon-number-resolving 
detection devices. One or both photons from a pair may 
be lost before they reach their detection devices. The 
reasons may be geometric filtering (one photon from a 
photon pair is not collimated to the detector area), re- 
flections on optical elements in experimental setup or ab- 
sorption of a photon along its path to a detection device. 
We describe this effect by two beamsplitters BSs and 
BS/ [9] placed in the signal-field and idler-field paths, 
respectively. We model a photon-number-resolving de- 
tection device as a multiport lxJV [10] followed by N 
single-photon detectors. This description is applicable 
also when an intensified CCD camera is used (as in our 
experiment). We note that detectors able to resolve di- 
rectly photon numbers to some extent have been con- 
structed [11-13]. From practical viewpoint, detectors us- 
ing time multiplexing (reached in fiber optics) and one 
or two single-photon detectors are promising [14-16]. We 
note that an intensified CCD camera has already oc- 
curred to be useful when studying spatial correlations 
of photon pairs [17]. 

We assume that the signal and idler fields in the out- 
put plane of nonlinear crystal NLC are described by the 
following statistical operator psi written in Fock basis: 



oc oo 

psi = E E P( n s, n i)\ n s)ss{ns\ ® (1) 

ns- rii— 

the symbol p(ns,ni) denotes the joint signal-idler 
photon-number distribution. 

Statistical operator pgj appropriate for the signal and 
idler fields in front of the detection devices can be written 
as: 

oo oo 

Psi = E E Pfas.ni) 

ns—0 ni=0 

x E (T s ) T s R s S ~ ls \ l s)ss(ls\ 

ls=0 ^ ' 

xE(";)t^ j R-^|Z / ) jj (/ / |. (2) 

The symbols Rs and Ri (Tg and Tj) denote intensity 
reflectivities (transmissivities) of the beamsplitters in the 
signal-field and idler-field paths. 

We assume a multiport 1 x Ns (1 x Nj) followed by 
Ns (Ni) single-photon detectors with quantum efficien- 
cies rjSi (riii ) and dark-count rates dsi (di i ) in the signal 
(idler) path. Detection of a photon in the fc-th detector 
is described by the detection operator Dk [6]: 

oo 

Dk = E^ 1 - (! - Vk) n ] + 4(1 - Vk) n }\n) kk (n\. (3) 

n=0 

Detection operator _D£° corresponds to the case when no 
detection has occurred: 

W = 1 - Dk- (4) 

The effect of splitting photons in the signal field in 
a 1 x Ns multiport can be described by the relation 
as = yi^—i tSi O'St ; where as is the annihilation operator 
of the signal field entering the multiport, the annihila- 
tion operator as i describes a field at the i-th multiport 
output, and ts t stands for amplitude transmissivity of a 
photon from the input to the i-ih output. 1 x JV; multi- 
port in the idler-field path is described similarly and the 
symbol tj i then refers to amplitude transmissivity of a 
photon from the input to the i-th multiport output. 

The probability Cgo ,/d that given cs detectors in the 
signal field and given c/ detectors in the idler field detect 
a photon whereas the rest of detectors does not register 
a photon is determined as follows: 

C s ojo = Tr s/ J J] D a J] DF 
{ aes D bes\s D 

x n d < n d t | • ( 5 ) 

cgi d dei\i D J 
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The symbol S (I) denotes the set of all signal-field (idler- 
field) detectors S = {Si, S Ns } (I = {h,..., InA)- 
The set S D (I D ) contains signal-field (idler-field) detec- 
tors that registered a photon. 

Using the statistical operator pgj in Eq. (2) the prob- 
ability C s d jd defined in Eq. (5) is written in the form: 
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having cs detections somewhere at N$ signal detectors 
and cj detections somewhere at Nj idler detectors can 
be expressed as: 



/ 
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Using the expression for C s ojd in Eq. (6) we arrive at 
the relation: 
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If the number of photons detected by an intensified 
CCD camera is much lower than the number of pixels 
detecting the field with a nonnegligible probability, the 
limits Ns — ► oo and Nj — ► oo are appropriate. When 
determining these limits, the overall noise levels Ds and 
Dj are hold to be constant (Ds = Nsds, Dj = Nidi). 
The coefficients K defined in Eq. (9) then simplify: 



min(ci,ni) , 

K i '°°{c i ,n i ) = E \l 

D Ci-l 
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(T^)'(l-T^)"'-' 
exp(-A), i = S,I. (10) 
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We now consider symmetric multiports (ts 1 = ts 2 = 
■ • • t SNs =ts,ti 1 =ti 2 = ... = ti Ni = ti) and detectors 
having the same characteristics in the signal and idler 
fields (r) Sl =rjs 2 = ... = vs Ns = Vs, d Sl = d S2 = . . . = 
ds Ns = d s , rjh = Vh = ■ ■ ■ = Vinj = Vi, d h = d i 2 = 
... = di Nj = di). Then the probability f Ns - Nl (c s , cj) of 



III. RECONSTRUCTION OF THE JOINT 
SIGNAL-IDLER PHOTON-NUMBER 
DISTRIBUTION 

The probabilities (frequencies) /°° ,00 (c5, Cj) are mea- 
sured and the relation in Eq. (8) then has to be inverted 
in order to obtain the joint signal-idler photon-number 
distribution p(ns,ni). The relation in Eq. (8) together 
with the coefficients K l '°°(ci,rn) defined in Eq. (10) can 
be inverted under special conditions analytically. For in- 
stance, if Ds = Di = the inversion relation is found 
using "convolution" of f°°<°° with the Bernoulli distribu- 
tions with efficiencies l/(Tst]s) and 1/(7/7//). However, 
analytical approaches are not suitable for processing real 
experimental data [18]. Reconstruction algorithms have 
occurred to be suitable in this case [14]. Such algo- 
rithms are able to find a joint signal-idler photon-number 
distribution p^°°\ns,ni) that gives the measured fre- 
quencies f°°'°°{cs,ci) with the highest probability from 
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all possible photon-number distributions. We use the 
Kullback-Lciblcr divergence as a measure of distance be- 
tween the experimental data and results provided by the 
developed theory. The joint signal-idler photon-number 
distribution p(°°\ns,ni) that minimizes the Kullback- 
Leibler divergence may be found, e.g., using the iterative 
Expectation-Maximization algorithm [19,20]: 

/ 00 ' °(is ) i/)^ s ' 00 (is,ns)^ 7 ' 00 (i/ ) n 7 ) 



(a) 



v x 

i s ,ii=0 Z-'isjl- 



:0 K s >°°(i s ,js)K I >°°(i I ,j I )p( n )(j s ,ji)' 

(11) 



The symbol p^ n \ns,ni) denotes the joint signal-idler 
photon-number distribution after the n-th step of iter- 
ation, p(°\ns,ni) is an arbitrary initial photon- number 
distribution. 

Experimental setup used for the measurement of fre- 
quencies f°°'°°(cs,ci) is shown in Fig. 2. 
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FIG. 2. Scheme of the setup for detection of photon pairs. 
Fields composed of typically tens of photon pairs are gener- 
ated in a nonlinear crystal, then they propagate through a fre- 
quency filter and are detected on an intensified single-photon 
CCD camera. 

Photon pairs are generated in a 5 mm long LHO3 crys- 
tal pumped by intense ultrashort pulses delivered by a 
titanium-sapphire femtosecond laser followed by a regen- 
erative amplifier at the wavelength of 800 nm. The laser 
system runs at the repetition rate of 11 kHz and after 
converting the 800-nm beam to its second harmonic it 
typically delivers pulses with the energy of 0.1 p,3. Pho- 
tons in the idler field are reflected from a high-reflectivity 
mirror and impinge on an intensified CCD camera. Sig- 
nal photons come directly to the photocathode of the 
intensified camera. Three regions of interest are defined 
in the field of view of the camera: the first one is for 
the signal field, the second one for the idler field, and 
the third region serves for reference measurements of the 
noise level. The whole field of the camera is filtered by a 
high-transmittance high-pass filter blocking light below 
700 nm and by an intereference filter of 10 nm FWHM 
centered at 800 nm. The interference filter selects nearly 
degenerate photon pairs. A typical histogram of the mea- 
sured frequencies of the joint signal-idler photon-number 
distribution f°° ,0 °(cs, cj) is shown in Fig. 3. 
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FIG. 3. (a) Measured histogram of the frequencies J 00 ' 00 as 
a function of the signal (ns) and idler (m) photon numbers, 
(b) Topological graph of the histogram /°°'°°. 

The reconstruction algorithm in Eq. (11) lead 
to the joint signal-idler photon-number distribution 
p(°°\ns, fii) in the output plane of the crystal as shown 
in Fig. 4. We have assumed Tgrjs = Tjr\i = 0.03 af- 
ter taking into account losses in the setup (frequency 
filters, reflection on the output plane of the crystal, fi- 
nal quantum efficiency of the intensified CCD camera). 
The measurement has also provided values of noises; 
Ds = Dj = 0.1. The initial joint signal-idler photon- 
number distribution p^ was assumed to be uniform. 
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FIG. 4. (a) Reconstructed joint signal-idler photon-number 
distribution p^°°^ as a function of the signal (ns) and idler (nj) 
photon numbers, (b) Topological graph of the distribution 
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The graph in Fig. 4 clearly shows that signal and idler 
photons are generated in pairs to some extent. In order 
to quantify this property we use covariance C p of two 
stochastic signals m and ni described by a joint proba- 
bility distribution p(n\, 712): 
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Covariance C/ of the measured distribution J 00 . 00 equals 
0.025 in our experiment. This value of Cf is in agreement 
with the theoretical value determined along the model de- 
veloped in Sec. 2 for values of parameters defined above. 
The obtained low value of Cf shows that photons from 
photon pairs are lost with a high probability owing to 
low detection efficiencies. However, covariance C p of the 



reconstructed joint signal-idler photon-number distribu- 
tion p(°°^ equals 0.81. This means, that the reconstruc- 
tion algorithm is able to reveal clearly the original pairing 
of photons. We note that we have checked that the re- 
construction algorithm cannot give correlations of signals 
being independent (i.e. with no correlation) in the input 
to this algorithm. Independent simulation data provided 
by the theoretical model in Sec. 2 were used for this pur- 
pose. There are two reasons why the value of C p char- 
acterizing the reconstructed photon-number distribution 
p(°°) docs not approach value 1. Impossibility of the de- 
veloped model to describe precisely all noises occurring 
in the experiment represents the first reason. The second 
reason lies in the fact that the numerical implementation 
of the reconstruction algorithm looses its precision with 
the increasing number of iterations. 

Pairing of the signal and idler photons results in nar- 
rowing of the distribution p^^ of the difference ns — Jij 
of the signal and idler photon numbers. Figure 5 shows 
the distribution p^ determined from the reconstructed 

p(°°^(ns, ni) and compares it with the distribution p mdcp 
characterizing independent marginal signal and idler dis- 
tributions: 
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The symbol 6 stands for Kronecker delta. 
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FIG. 5. Distributions p_ (solid curve) and p™ dep (dashed 
curve) of the difference ns — ni of the signal and idler photon 
numbers. 
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We note that narrowing of the distribution p^°^ has 
been observed for stronger fields generated by parametric 
downconversion in a resonator utilizing correlations of 
photocurrents from two detectors [21]. 

Narrowing of the distribution pf ^ is accompanied by 
broadening of the distribution p+°^ for the sum ns+ni: 



(°°) / \ 
P+ » 



E E 

ns—O ni—0 



p (00, (n s ,n/). (15) 



This stemms from the following expression for the vari- 
ances ([A(ns - ni)} 2 ) and ([A(ns + ni)] 2 ): 

([A(n s ±n I )] 2 ) = ({An s ) 2 ) + ({An I ) 2 ) 

± 2^((An s ) 2 )((Anj) 2 )C p ; (16) 

covariance C p is given in Eq. (12). Figure 6 shows the 
distribution p^ 00 ' derived from the reconstructed data to- 
gether with the distribution p^ dcp determined from the 



independent distributions p^ and pf°' (14): 
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FIG. 6. Distributions (solid curve) and p^ dcp (dashed 
curve) of the sum ns + ni of the signal and idler photon 
numbers. 

Statistics of generated photon pairs as they occur in 
the output plane of the crystal cannot be inferred di- 
rectly from values of diagonal elements of the recon- 
structed joint signal-idler photon-number distribution 
p(°°) i n Eq. (ii) because of its blurring. This statistics 
can be determined from the marginal photon-number dis- 
tributions /Og°°' ) and or from the photon-number dis- 
tribution of the mixed signal and idler field p+°^ . Statis- 
tics are then judged according to the value of the coeffi- 
cient S, 



bp ~ (n) 2 (n) ' 



(18) 



i k ) = J2n k p(n),k = 1,2; 



S p — characterizes Poissonian statistics whereas S p = 1 
stands for Gaussian statistics. The values Sf.s — 1.001, 
Sfj — 1.007, and Sf,+ = 1.019 characterize the directly 
measured distribution fi 00 ' 00 ). The reconstructed distri- 
bution 

p(°o) i s 

characterized by the values S Pt s — 1.003, 
S p j = 1.013, and S Pt+ = 1.027. 

Assuming Poissonian statistics of photon pairs, i.e. 
p(ns,ni) = 5 ns ^ ni p ns exp(— p)/n s l, the theoretical 
model gives for p = 20 S p .s = 1, S p j = 1, S Pi+ = 1.025 
in the output plane of the crystal, whereas Sf y s = F 
Sfj — 1, = 1.017 are appropriate for distribu- 

tions detected by an intensified CCD camera. On the 
other hand, Gaussian photon-pair statistics, p(ns,ni) = 
<Sns,n7M ns /(M+l)" s+1 ! use d in the theoretical model with 
p. = 20, give Sp.s — 2, S p j — 2, S p .+ = 2 in the out- 
put plane of the crystal and Sf t s = F69, Sfj = 1.69, 
Sf,+ = 1.71 for distributions detected by an intensified 
CCD camera. Comparison of these theoretical results 
with the above ones stemming from experiment and re- 
construction clearly shows that photon pairs are gener- 
ated with Poissonian statistics. 

Poissonian statistics correspond to the physical pic- 
ture in which photon pairs can be generated into a great 
number of independent modes distinguishable in space 
and time (see Appendix A). In case of our experimen- 
tally generated fields containing typically several tens of 
photon pairs, each pair occupies its own mode with a 
high probability. If photon pairs are generated into one 
well-defined space-time mode, statistics of photon pairs 
would be Gaussian (see Appendix A). 



IV. JOINT SIGNAL-IDLER 
INTEGRATED-INTENSITY DISTRIBUTION 

A complete characterisation of the generated signal 
and idler fields is given by the joint distribution of their 
integrated intensities Ws and Wj in the output plane of 
the crystal; 



W, 



-f 



drEt\r)Er(r), 



(+)/ 



l = SJ. 



(19) 



The operator amplitudes E^ and e\~^ in Eq. (19) are 
given by Eq. (A5) in Appendix A. Integration over time 
variable in Eq. (19) is carried out over the field generated 
by one pump pulse. 

Inverting the photodetection equation the joint signal- 
idler integrated-intensity distribution P(Ws,Wi,s) re- 
lated to s— ordering of field operators can be determined 
using the joint signal-idler photon-number distribution 



as ( [22], chap. 4): 
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L° are Laguerre polynomials. The parameter s equals -1, 
0, and 1 for antinormal, symmetric, and normal ordering 
of field operators, respectively. The intcgratcd-intensity 
distribution related to normal ordering does not include 
noise from vacuum fluctuations and that is why it shows 
all features of the generated fields in their full complex- 
ity. However, its structure is rather complex in our case 
of nonclassical correlated fields (generalized functions are 
necessary for the description) and that is why we use the 
distribution related to symmetric ordering. This distri- 
bution behaves well and still has the capability to resolve 
nonclassical properties of the fields by its negative values. 
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FIG. 7. Joint signal-idler integrated-intensity distribution 
P(Ws, Wi, 0), (a) 3D plot, (b) topological plot showing lines 
having positive values, and (c) topological plot showing lines 
having negative values of the integrated intensity for the re- 
constructed photon-number distribution p'^^ng, n/). 

The signal-idler integrated-intensity distribution 
P(Ws, Wi, 0) related to symmetric ordering of operators 
and obtained from the reconstructed signal-idler photon- 
number distribution shown in Fig. 7. It clearly 
shows strong correlations in values of the integrated in- 
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tensities Ws and Wj. This correlation inherent to the 
fields in the moment of their generation again shows that 
photons are generated in pairs, i.e., the same number of 
energy quanta is put into both fields during their gener- 
ation. Negative values of the distribution P(Ws,Wi,0) 
(see Fig. 7c) reached in some regions of values of inte- 
grated intensities Ws and Wi mean that the generated 
signal- idler field with its correlations is nonclassical (i.e., 
it cannot be described by classical statistical optics) . For 
comparison, Fig. 8 contains the signal-idler integrated- 
intensity distribution P(Ws, Wj, 0) for an ideally corre- 
lated signal- idler field with Poissonian statistics. 



(a) 




FIG. 8. Joint signal-idler integrated-intensity distribu- 
tion P(W s ,Wi,0), (a) 3D plot, (b) topological plot 
showing lines having positive values, and (c) topo- 
logical plot showing lines having negative values for 
p^\n s , ni ) = <5„ s , n/ M" s exp(- M )/n s !, H = 20. 

All nonclassical features visible in this distribution in 
Fig. 8 are contained in the experimentally-obtained dis- 
tribution shown in Fig. 7. 
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V. CONCLUSIONS 

We have experimentally verified that photons are gen- 
erated in pairs into the signal and idler fields in the non- 
linear process of spontaneous parametric downconver- 
sion. Both the experimentally determined joint signal- 
idler photon-number distribution and joint signal-idler 
integrated-intensity distribution show strong correlations 
between quantities characterizing the signal and idler 
fields. For example, covariance of the signal and idler 
photon numbers is greater than 80 %. Statistics of the 
generated photon pairs have been determined to be Pois- 
sonian. Nonclassical character of the generated signal- 
idler field manifests itself by negative values of the joint 
signal-idler integrated-intensity distribution. 
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APPENDIX A: PHOTON-PAIR STATISTICS IN 
SPONTANEOUS PARAMETRIC FREQUENCY 
DOWNCONVERSION FROM FIRST 
PRINCIPLES 

The interaction Hamiltonian fl in t of the process of 
spontaneous parametric frequency downconversion can 
be written in the form [1,23,24]: 

H int (t) = ant 12 12 ^ k s> *)4(m4(M + h.c. 

= + (Al) 

where C; n t is an interaction constant and <fr(ks, fcj, t) de- 
scribes a detailed space-time structure of the emitted 
photon fields. The symbol H.c. denotes Hermitian con- 
jugate. The operator (H^) stands for the part 
of the interaction Hamiltonian fl; n t containing creation 
(annihilation) operators of modes in the signal and idler 
fields. 

The state of the signal and idler fields in the output 
plane of the crystal is obtained solving the Schrodinger 
equation in the form: 



71 = 



|V-o) = |vac), 



|V>n) 



Tl p OO rT\ 
oo J —OO 



X flint (Tl) . . .fli„t (T n ) | Vac), 



dr n 



n= 1,2 



(A2) 



The signal and idler fields are assumed to be in the vac- 
uum state | vac) in the input plane of the crystal. 

If the number of photons in the signal and idler fields 
is much lower than the number of independent modes 
constituting these fields, we may approximately write: 

'^bU L dn L dT2 -L dTn 



fl^Vi)---^V«)|vac) 



int 

i 

h 



| vac), 



1,2, 



(A3) 



The state \i() n ) then describes the field containing just n 
photon pairs. 

Statistics of the generated photon pairs may be in- 
ferred from the normally ordered moments of the "cre- 
ation operator of photon pairs" fl pa i r and its Hermitian 
conjugate: 



P pair ( Ts , tj) = 4 \t s )E\ >(tj) . 

?(+) 



*(-)/ 



(A4) 



The symbol E^' (e[ ^) denotes the positive- (negative-) 
frequency part of the clcctric-ficld amplitude of the signal 
(I = S) and idler (I = I) field: 

4 (+) ( r ) = ^2e l (k l )ai(k l )c^p(-iuj kl T), (A5) 



E, 



(-) 



-( 



ei(ki) is a normalization amplitude of the mode ki, I = 
S, I. Underlining of the operators on the right-hand side 
of Eq. ( A4) means that "only the signal and idler photons 
created in the same elementary event are considered" (see 
the expression for \ip n ) in Eq. (A3)). 

We may write for the normally ordered moments of 
^pak an d flpair in the framework of the above defined 
approximation: 



3=1 



flpair (T Sn+i-j , T I n + 1 _ j ) 

3=1 



II^I^airK.^JIvac) 
3=1 



n ^ va ° i ^ air ( ts -+ i -v i -i ) i ^i ) 



3=1 

nl^l^airK.r/Jlvac) 2 , 

3=1 

n=l,2,.... 



(A6) 
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Assuming t Si = ts 2 = . . . = t s and t Ti = r/ 2 = . . . = rj, 
Eq. (A6) provides a simplified expression: 



^pair (TS,T/) \lp n ) = 



(HPLA T s,Ti)\v&c) 



In 



(A7) 



Assuming that the contribution from the state \tp n ) is 
much greater that those from the states \ipk) for k = 
n + 1, n + 2, . . ., the relation in Eq. (A7) implies that 
statistics of generated photon pairs is Poissonian. 

If photon pairs are generated into one well-defined 
space-time mode, than the positive-frequency part of the 
electric-field amplitude of field I can be written as 



e[ +) (t) =a l (r)a l , 



(A8) 



where a/ describes amplitude per photon and di stands 
for annihilation operator of a photon in field I (I = S,I). 
We then have: 



(TS j Tl) 



|a(r s )a(r / )| 2 "(V„| 



44 



pair 

[fflSflf 1 " \lpn) ; 



|a(r s )a(r / )| 2n n! = n! (V-il^fo, Tr)|vac 



2n 



(A9) 



Relations in Eq. (A9) are in agreement with Gaussian 
statistics of generated photon pairs. 
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